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techniques are poorly suited for analyses of sources with extensions comparable to the detector’s 
field of view (greater than ~1°). Nearby pulsar wind nebulae, supernova remnants interacting 
with molecular clouds, and dark matter signatures from galaxy clusters are just a few potentially 
highly extended source classes. A three dimensional maximum likelihood analysis is in devel¬ 
opment that seeks to resolve this issue for data from the VERITAS telescopes. The technique 
incorporates relevant instrument response functions to model the distribution of detected gamma- 
ray like events in two spatial dimensions. Additionally, we incorporate a third dimension based 
on a gamma-hadron discriminating parameter. The inclusion of this third dimension significantly 
improves the sensitivity of the technique to highly extended sources. We present this promising 
technique as well as systematic studies demonstrating its potential for revealing sources of large 
extent in VERITAS data. 
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1. Introduction 

Imaging atmospheric Cherenkov telescopes (lACTs) have had great success in detecting sources 
at TeV energies. However, the sources typically detected are point sources or exhibit only mod¬ 
erate extensions (radial extension < 0.5°). Other instruments with larger fields of view have de¬ 
tected sources with very large extensions (>1°) including the cocoon detected by Fermi-LPIV in 
the Cygnus region, believed to be the result of freshly accelerated cosmic-rays confined wifhin a 
cavify in fhe Cygnus X sfar-forming region [1]. The morphology of fhe cocoon was characferized 
using a Gaussian wifh a = 2.0° ± 0.2 making if highly extended. Milagro also reported evidence 
of exfended emission (infrinsic full widfh af half maximum of 2.6°1 q g) [2] coincidenf wifh fhe 
Thrmf-LAT Geminga pulsar [3]. This extended source, idenlified as MGRO J0632-I-17, is fhus po- 
fenfially a pulsar wind nebula (PWN). Addifionally, dark mailer signalures from galaxy clusters can 
exlend over regions larger lhan a degree [4]. Detecting Ihese sources requires an analysis which 
can accounf for all of fhe emission in fhe region simullaneously. This is nol a Irivial lask given 
fhe analysis techniques used in lACT dafa analysis. Sfandard melhods for eslimaling background 
emission for Ihese inslrumenls are eilher fhe ring background model (RBM) or fhe reflecled-region 
model [5]. 

The RBM analysis estimates fhe level of background for a particular region of inleresl {ON 
region) by using an annulus, or ring, around fhe region. In order lo obfain an accurale measure of 
fhe background fhe ralio of fhe ON region lo fhe ring region is typically on fhe order of ~l/7. The 
reflecled-region model requires observalions lo be faken offsel from fhe source posilion. Suilable 
OFF regions are chosen af posifions in fhe field of view which have similar offsels from fhe fracking 
position as fhe eslablished ON region. For sources wifh exfensions >0.5°, no OFF regions can be 
determined and a larger observing offsel is needed, reducing sensilivily fo fhe source. Ulfimafely, 
fhe field of view of fhe inslrumenf limils fhe size of source which can be observed wifh fhis tech¬ 
nique. In fhe RBM highly exfended sources may require a larger ON fhan OFF region, resulting 
in a poorly defermined background. These examples for fhe RBM and reflecled-region analyses of 
sources wifh large exfensions are predicafed on fhe assumplion lhaf fhe morphology of fhe source is 
known. For fhe mosf exfended sources fhe morphology is rarely undersfood wifh absolute cerfainfy. 
In fhis case, defermining suilable OFF regions becomes very difficulf and a cerfain amounf of self 
sublraclion due lo source contamination in the OFF region is almost guaranteed. Thus, a different 
analysis technique is necessary to detect these sources in VERITAS data. 


2. VERITAS 

The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is an array of four 
12-meter lACTs located at the Ered Eawrence Whipple Observatory (EEWO) in southern Arizona 
(31°40’ N, 110°57’ W, 1.3km a.s.l.). The instrument covers an energy range from 85 GeV to > 30 
TeV with an energy resolution between 15 - 25%. VERITAS is capable of detecting a point source 
with 1% the Crab nebula flux in ~25 hours and has an angular resolution (68% containment radius) 
of <0.1° at 1 TeV. 
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3. VERITAS Maximum Likelihood Method 

In order to tackle extended sources in VERITAS data, a three dimensional maximum likeli¬ 
hood method (MEM) is in development. In general, the method of maximum likelihood works by 
modeling the expected distribution of events in data. A likelihood value L can then be computed 
from these models according to 

L(s)=np(e,-|s) (3.1) 

i=\ 

where p is the model, i indexes individual data events, 6 is a set of observables which represent the 
data, and s is a set of free parameters chosen to maximize L. Typically, the log of the likelihood is 
computed as this turns the product into a sum and the computation is simplified. 

The VERITAS MEM models consist of both a two dimensional spatial model and a mean 
scaled width model (MSW) for both the y-ray (source) and background components. In this way 
the full log-likelihood for a single observation becomes 

d 

~ ^obs^OgfNexp) Nexp~\~ [Sirc(r/|s) ' W(src(lVj|s) -|- Shkgdifi) ' (3-2) 

/=! 

To correlate the spatial source model {Ssrc) with the MSW source model (1^^^). the product of the 
two distributions is taken (similarly for the background component). An additional term is included 
in the likelihood to model the total expected number of events from both signal and background. 
The expected number of events in a dataset is taken to be Poisson distributed with expected events 
Nexp and observed events Nobs- The first two terms in Eq. 3.2 are based on this assumption. The 
data is treated as unbinned in the spatial and MSW dimensions and coarsely binned in reconstructed 
energy. The free parameters are spectral parameters associated with the source being modeled. 

It follows that to analyze data from multiple observing positions, each field can be indepen¬ 
dently modeled and the resulting log-likelihoods summed to get a likelihood for the entire observa¬ 
tion set simultaneously. Since all models share the same free parameters (which are independent of 
the observing conditions between fields) this can be extended to observations which cover a range 
of observing positions, detector configurations, or even to connect observations from different de¬ 
tectors. 


4. Defining the Input Models 

In this section, the various models which represent the expected distribution of the data, and 
how they are derived is detailed. Specific attention will be paid to the source spatial model in 
Section 4.2 outlining the details of how we incorporate the various instrument response functions 
(IREs) from VERITAS. Each model component is derived from either data (primarily for modeling 
background) or simulated y-ray air-showers (for source models). In both cases, events were pro¬ 
cessed using standard medium cuts for all variables with the exception of MSW in which a larger 
range was used to better characterize the background models. 
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4.1 Background Spatial Models 

The background spatial models represent the expected distribution of y-ray-like events across 
the field of view. They are derived from observations taken on weak blazars and dwarf spheroidal 
galaxies. Data taken in the galactic plane is not used to prevent possible contamination from diffuse 
y-rays. These data are first transformed into the spatial coordinate system used in the analysis and 
then binned into 0.01° spatial bins. A radial acceptance distribution is generated in by summing 
all bins within a range of offsets from the field of view center. This radial acceptance is extrapolated 
into two dimensions to form the spatial model associated with the expected background. 

Because data is used to derive the background spatial models there are two main effects which 
must be accounted for. The first effect is contamination from any potential source within the field 
of view which results in an overestimate of the background at the offset of the source. Secondly, 
stars in the field of view can cause deficits in the surrounding region. To alleviate these effects all 
bins within a set angular distance* of known (or potential) sources and bright stars^ are excluded. 
To correct for the excluded regions, each bin of the radial acceptance curve is weighted by the 
relative number of bins which actually contribute to it. 


4.2 Source Spatial Model 


In order to produce a y-ray source spatial model the formulation from Mattox et al. (1996) Eq. 
2 is used [6]. 


5^rc(f|s) 




* E(r , E')]S{E'\s)R(f, E',E)A (r, E')dE'dE 


(4.1) 


This gives the expected distribution of the source emission as it is observed. B represents the intrin¬ 
sic source morphology, which in the case of a point source is a delta function centered at the source 
position. The various instrument response functions (IRFs) accounted for are the point spread func¬ 
tion (PSF, P), energy dispersion (/?), and effective collection area (A). Additional parameters are the 
intrinsic source energy spectrum {S), the spatial coordinates in the field of view (r), the true y-ray 
energy (£"), and the observed y-ray energy {E). To make the computation of this model tractable, 
certain assumptions are employed. These include finely binning B and P spatially as well as finely 
binning in The spatial model then becomes a summation rather than an integration 


“^ircCrij'Is) — Bffi fiPffi fi(Ti j')Afji fi(^Ej^') I Rffifi[E,Ej^')dE 

k,m,n •'Emin 


nF' 

I kdipper 


F' 

‘^k.lower 


S[E'\l)dE'. 


(4.2) 


This now gives the value of the spatial model within a given bin i,j of the spatial map. Potential 
“smearing” of events due to the instrument resolution is accounted for by summing the contri¬ 
butions from all other bins n,m at position i,j in the spatial map. A pictoral example of the input 
source model and resulting MFM computed source model for both a point source and 1.5° diameter 
extended disc source are given in Figures 1 and 2 respectively. 

This formulation solves an additional challenge. Following Eq. 4.1 the entire source model 
and all relavent contributions from each IRF need to be recomputed at each step of the fit. This 

* The size of the exclusion region for models in this paper is 0.4°, however this value is still undergoing optimization. 
^bright stars are defined as those with B magnitude 8 or helow. 

and P are binned spatially using 0.025° x0.025° bins. E’ is binned in logioiE’ [TeV]) at increments of 0.05. 
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(a) Input point source model (b) Computed point source model 

Figure 1: Example of a computed point source model using Eq. 4.2. (a) Input source model as¬ 
suming a delta function at the source position (~0.5° offset from the camera center), (b) Resulting 
source model as computed by the MEM described here. Eor this test case, a Crab-like source 
spectrum was used. 



(a) Input extended disc model (b) Computed extended disc model 

Eigure 2: Example of a computed extended source model using Eq. 4.2. (a) Input source model 
assuming a 1.5° diameter disc centered at 0.5° offset from the camera center, (b) Resulting source 
model as computed by the MEM. Eor this test case, a 10% Crab-like source spectrum was used. 


is computationally intensive. However, using Eq. 4.2 the component of the source model derived 
from the IREs (which are independent of the source spectral parameters) can be precomputed for 
every bin of E’ requiring only a reweighting of these bins for each iteration of the fit. 

The IREs are parameterized for a range of observing conditions including zenith angle, az¬ 
imuth, night sky noise level, observing offset from the source, atmospheric conditions, telescope 
array configuration, and number of telescopes contributing to the shower reconstruction. The PSE 
is modeled based on fits of a King function to simulated y-ray air-shower events generated at fixed 
energies of 0.1, 0.3, 1, 3, and 10 TeV across fhe full range of observing condifions lisfed above. 

+ ' (43) 

In particular X and a are compufed for each combinafion of possible observing condifions and 
inferpolafed on fo obfain a PSE model af any pofenfial observing posifion and energy. These sim- 
ulafions are also used fo generafe fhe energy dispersion of fhe insfrumenf. As fhe inferpolafed 


5 































A Novel Method for Detecting Extended Sources with VERITAS 


Joshua V. Cardenzana 



Figure 3: Energy dispersion as a function of logarithmic true energy. The plot represents the 
fraction of y-ray events reconstructed at a given energy (Erec) as a function of (Erec-Etrue)/Efrue- 
Values are based on fixed-energy y-ray air-shower simulations at 70° elevation, 180° azimuth, and 
offset by 0.25° from the camera center. Simulations are processed using non-standard cuts on MSW 
and including all events within the field of view. Values are interpolated at energies not represented 
by the fixed-energy simulations. 

energy dispersion becomes less well defined below ~300 GeV (due to limited statistics and the 
interpolation method used, see Figure 3) the MEM is currently restricted to analysis above this 
energy. 

4.3 Mean Scaled Width 

In order to improve sensitivity to extended sources a third dimension based on a y/hadron 
discriminating parameter known as mean scaled width (MSW) [7] is incorporated. This parameter 
is based on the average of the Hillas width parameter [8] from images of all participating telescopes 
normalized by the expected width derived from simulated y-ray air-showers. 



(4.4) 


n is the number of participating telescopes and j represents the various shower parameters over 


which the expected width value is determined, such as integrated image intensity and impact dis¬ 


tance. It follows that y-ray events will have a MSW distribution which peaks very close to 1. 
However, cosmic-rays, which typically result in less compact images, will peak at higher values of 
MSW (see Figure 4). For the MEM analysis, MSW is restricted to the range 0.8 - 1.3 to provide 
ample distinction between the source and background models. 

The source MSW model is derived from simulations of y-ray air-showers while the back¬ 
ground models are derived from actual observations, as mentioned above^. 

5. Validation Plans 

The planned validation of the VERITAS MEM analysis will initially involve a self consistency 
check in which the generated models are used to produce a sample of toy Monte Carlo events. 

^ As there is potentially a correlation between MSW and spatial position, both sets of MSW distributions are derived 
for a range of different offsets. 
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Mean Scaled Width 


(a) Gamma-ray MSW Distribution (b) Background MSW Distribution 

Figure 4: MSW distributions derived (a) for the source model from simulated 7 -ray air-showers 
and (b) for the background model from observations. Distributions are shown as a function of 
logarithmic energy. The color scale represents the number of events used to generate the distri¬ 
bution. Models are extracted for an observation position of 70° elevation, 180° azimuth, and 0.5° 
offset with four telescopes participating in the shower reconstruction. The slight peak near 1 in the 
background MSW distribution below ~300 GeV is currently under investigation. 


These generated events will then be fit with the same models to ensure the fitted result is consistent 
with inputs. The next step will be to fit the models to dedicated OFF observations to check the 
background models are consistent with data. Once the background models have been verified, 
a test on the Crab Nebula, a known bright point source, will be conducted to check the source 
model generation is accurate. Additionally, data has been taken on the Crab Nebula following a 
raster scan pattern to produce an extended source [9]. Tests on these datasets will provide a check 
of the extended source fitting using a source with a known spectrum. Finally, the technique will 
be applied to known extended sources such as supernova remnant G106.3-I-2.7 [10] and MGRO 
J1908-1-06 [11]. Each check is designed to ensure all components of the MEM are understood and 
behave as expected as well as to identify regions where the technique struggles. 

6. Discussion 

The expected number of 7 -ray events predicted by the MEM for the Crab-like source model 
presented in Eigure 1 was compared to a sample of observations of the Crab nebula taken under 
similar observing conditions and analyzed with a standard RBM analysis. The MEM prediction 
was found to be consistent to within ~10%. Additionally, the extended source (which uses a 10% 
Crab-like spectrum) actually predicts ~8-9% of the expected Crab flux. This is due to the reduced 
instrument sensitivity at larger offsets from the center of the field of view which is reflected in the 
asymmetry of the computed model (see Eigure 2b). 

The analysis is still under development and will undergo a set of optimizations in the fu¬ 
ture. It will be tested whether extending the range on MSW provides better characterization of 
the background and thus improved separation between source and background model components. 
Additionally, how the chosen coarse binning in observed energy affects the result will be investi- 
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gated. Improvements to the VERITAS shower and energy reconstruction methods will also lead to 
more accurate IRE parameterizations and thus improve the overall MEM. 

Once validated, the technique can be employed on several known extended sources currently 
undetected by lACTs, such as the cocoon in Cygnus and MGRO J0632+17. This technique also 
allows the possibility for doing energy-dependent morphology fitting across multiple instruments 
using a set of common spatial templates. In so doing the VERITAS MEM could be an invaluable 
tool for multi-wavelength studies in the future. 
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